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1 . SCOPE 


This specification establishes the modifications to the CLASSY 
program as specified in IDSD Category 1 Job Order Task Agreement , 
titled CLASSY Program Modification. This modification includes 
the CLASSY post-processor MAXLABEL and its associated program 
IMAGE. 
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2 . APPLICABLE DOCUMENTS 


The following documents form a part of this specification: 

J. 0, 76-662 Task Agreement titled: CLASY Program Modification 

Technical Memorandum Mathematical Description and Program 
Documentation for CLASSY, An Adaptive Maximum Likelihood 
Clustering Method, by R. K. Lennington and M.E. Rassbach 
Elogic, Inc. Houston, Texas. 

Earth Observations Division Version of the Laboratory for 
Applications of Remote Sensing System ( EOD-LARS YS ) User Guide 
for the IBM 370/148 Volumn II - User's Reference Manual. 


3. SYSTEM DESCRIPTION 


The CLASSY clustering algorithm serves to estimate the component 
distributions which make up the overall mixture distribution of 
the data. The present software uses the statistics for these 
component distributions (clusters) along with a set of labelled 
data vectors to produce maximum likelihood estimates of the 
proportion of each labelled class associated with each cluster. 
These estimates (called betas) are obtained through a fixed 
point iteration procedure. The estimated betas are used to 
obtain two different estimates of the proportion of each labelled 
class in the whole scene. In addition labelled cluster maps 
are produced for the whole scene using two different techniques. 

The maximum Likelihood System of programs consists of the following 
four programs: CLASSY, DOT^ IMAGE and LABEL. (Figure 3-1). 

The CLASSY Program creates statistically meaningful clusters 
from the pixels in a selected area of the segment image tape. 

The program options are defined in a control card file. The 
cluster statistical parameters for each iteration are saved on 
the CLASSY Cluster Statistical Parameter Files for later use 
by the LABEL program and a CLASSY One Channel Unlabelled Cluster 
Map is optionally written to tape for use by the Accuracy 
Assessment Programs. Cluster maps are written at the end of the 
last iteration and optionally written after each iteration 
(Figures 3-1 and 3-2) . 

The DOT program creates a Ground Truth Selected Pixel Values and 
Names File associated with pixel radiance values from the same 
segment image tape used by CLASSY, the Ground Truth Description 
File and the program options as defined on the Dot Control Card 
File. (Figures 3-1 and 3-5) . 


The IMAGE program creates a Selected Pixel File from the same 
segment image tape used by CLASSY and DOT. The LABEL Control 
Card File is used to define the options. (Figures 3-1 and 3-9) . 

The LABEL program processes the CLASSY Cluster Statistical 
Parameter File frojn CLASSY, the Ground Truth Selected Pixel Values 
and Name File from DOT and the Selected Pixel File from jlMAGE 
to create a One Channel Labelled Cluster Map on tape for use 
by the Accuracy Assessment Programs and the six following 
reports : 

Iteration Report 
Pixel Labelled Dot Map 
Unlabelled Dot Map 

Entropy and Posterior Probability for each Cluster Report 
Pixel Labelled Cluster Map 
Cluster Labelled Map 
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Processing Order - EXECS: CLASSY or CLASSVN PROGAAK: CLASSY 



SELECTED 
PIXEL FILE 










C?*«*V #>{»C 


GLODAL TXTL1H FORTRAN CMSLI8 

nr gsTORAtsc of: :mmo isur 2 

ftTVPE I»rF STUN 2* REQUIRED 

•EXIT 

IIP ft INDEX EO 2 ISK1P 31 
ftIF S3 EG 1 ft SKIP # 
ftSK *12-1 . . _ . 

TAPE FSF ftSK 
S CONTROL OFF 


»»*AJ r'AGE IS 
OF POOR QUALITY 


CLASSY EXEC FILE 


GETpISK JSC770 191 350 R 8/A PASS AUCOIN 
ft SEEP • 1 

-GET GETDISK TEMP 3M CLEAR DETACH 

ftlF IRETCUllE EQ 0 *UKIP 7 

ft IF ftDEEP LE 7 ftSKIP 3 ' 

ftTYPE NO 25 CYL TEMP DISK AVAILABLE 

CP 0 T 

ftEXIT 

CP SLEEP 5 MIN 
ft DEEP • ft DEEP ♦ I 
ft GOTO -GET 
0 DISK D 


GETDISK 

LARSYS 











FILEDEF 

FTO2F001 

DISK 

FILE FT02F001 

Al (LRECL 320 PLKSIZE 

320 

*‘ERM 

FILEDEF 

FT02F002 

DISK 

FILE F702F002 

A1 (LRECL 320 PLKSIZE 

320 1 

PERM 

FILEDEF 

3 TERM (PERM 










FILEDEF 

FT04F001 

DISK 

FILE FT04F001 

A4 (LRECL 3060 PLKSIZE 

3060 

PERM 

FILEDEF 

FT04F002 

DISK 

FILE FT04F002 

A4 (LRECL 3060 BLKSIZE 

3060 

PERM 

FILEDEF 

FT04F003 

DISK 

FILE FT04F003 

A4 (LRECL 3060 PLKSIZE 

3060 

PERM 

FILEDEF 

6 DISK 0UT6 LISTING V 

(PERM 








FILEDEF 

FT09F001 

DISK 

FILE FT09F001 

Al (LRECL 320 PLKSIZE 

320 PERM 

FILEDEF 

FTUF001 

TAF’l 

(BLKSIZE 

3060 

RECFM 

U 

den 

1600 

PERM 



FILEDEF 

FT11F002 

TAPI 

(BLKSIZC 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F003 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F004 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

FERM 



FILEDEF 

FT11F005 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDCF 

FTUF006 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 



FILEDEF 

FT11F007 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FTUF0P8 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F009 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F010 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDF.F 

FT11F01 1 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F012 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FTHF013 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F014 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F015 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FTUF016 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F017 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEr 

FT11F018 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FTUF019 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FrilF020 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FTHF021 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F022 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F023 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FTUF024 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F025 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

rTHF026 

TAPI 

(PLKSIZE 

2060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F027 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FTUF028 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F029 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT11F030 

TAPI 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 



FILEDEF 

FT12F001 

TAP2 

(PLKSIZE 

3060 

RECFM 

u 

DEN 

800 PERM 



FILEDEF 

FT12F002 

TAP2 

(PLKSIZE 

12356 RECFM U DEN 800 

PERM 



FILEDEF 

FT13F001 

DISK 

FILE FT13F001 

A4 (LRECL 320 PLKSIZE ; 

320 PERM 

FILEDEF 

FT14F001 

DISK 

FILE TTMFOOl 

PI (LRECL 3060 PLKSIZE 

3060 

PERM 


FIL FT16F001 DISK FILE FT16F001 A1 (LRECL 3060 RECFM U PLKSIZE 3060 PERM 
F1LEDEF FT17FOOI TAP3 (BLKSIZE 12356 RECFM U HEN 600 PERM 


FILEDEF FT19F00I DISK FILE F119F001 Al (LRECL 1660 PLKSIZE 1860 PERM 

FILEDEF FT19F002 PISK FILE FT19F002 Al (LRECL I860 PLKSIZE 1660 PERM 

FILEDEF FT20F001 DISK FILE FT20F001 A4 (LRECL 320 PLKSIZE 320 PERM 

FILEDEF FT20F002 DISK FILE FT20F002 A4 (LRECL 320 PLKSIZE 320 PERM 

FILEDEF FT21F001 DISK ftl CC A1 (LRECL 80 PLKSIZE BOO PERM 
FILEDEF FT21F002 DISK FILE FT21F002 Al (LRECL 80 PLKSIZE 800 PERM 
FILEDEF 22 DISK FILE FT22F001 P (LRECL BOO PLKSIZE BOO XI ENT 2100 PERM 


FIL FT23F00! DISK SI 1 A4 (LRECL 120 PLKSIZE 130 RECFM W3 
FIL FT23F002 DISK Si 2 A4 (LRECL 120 PLKSIZE 130 RECFM OS 


FIL FT23F003 DISK ftl 
FIL FT23F004 DISK II 


3 A4 (LRECL 120 PLKSIZE 130 Rf.CFM VS 

4 A4 (LRECL 120 PLKSIZE 130 RfCFM VS 

FIL FT23F005 DISK ftl 5 A4 (LRECL 120 PLKSIZE 130 RECFM VS 

FILEDEF 26 DISK FILE MEANS A (PERM 

LOAD CLASY PLKCLA (CLEAR START NOMAP 

REMOTE E TO HOUSTON 

SF'OOL E HOLD 

PRINT OUT 6 LISTING D 

IIP ft INDEX LT 4 I SKIP 10 

CP DET 1U1 

TAPMOUNT ft 4 TAPI RI 

ft IF ft INDEX NE 5 ISK1P 3 

SIF 45 EG 1 ft SKIP 2 

ftSP -15-1 

TAPE F5F ISP 

FIL INMUVC DISK FILE FT16F001 A (PLKSIZE 3060 RECFM U FERM 
FIL OUT MOVE TAPI ( PLKSIZE 3060 RfcXf M U DEN 1600 PERM 
MOVEFILE 
I API. W1M 2 


PERM 

PERM 

PERM 

PERM 

PERM 


r;r him m 
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SAMPLE CLASSY CONTROL CARD FILL 


>TYPE 1394 CC A 


COMMENT 

CLASSY RUN ON SEGMENT 1394 

CHANNELS 

1 f 2 f 3 f 4 f *li f A f 7 f S 

DATAFILE 

FILE=1 

ITER 

5 

MAP 

5 

LINES 

1-11 7 

NPTS 

o 

Am 

SEGM 

1394 

PROC 

30f If SO 

TIME 

ISO 

♦END* 

(l»l)f (lfl)f ( 196f 1 > f ( 19Ar 1 17 ) r 

♦END* 



(If 117) 


>d*f stop 2m 

STORAGE * 02048K 

Rf T-O.Ol/O.Ql 10*45 • 37 


SAMPLE CLASSY EXECUTION of Segment 1394 


>eUssy 1394 5123 1 

GLOBAL TXTLIB FORTRAN CMSLIB 

JSC770 191 HAS BEEN ATTACHED AS 350* 

B (350) R/0 

350 HAS BEEN LOGGED IN AS B/A DISK * 

DASD 192 DETACHED 

YOU ALREADY HAD A VIRTUAL DEVICE 192. IT IS BEING DETACHED. 

TEMP 157 HAS BEEN ATTACHED AS 192. (003.00 MEGABYTES) 

192 HAS BEEN LOGGED IN AS D DISK. 

D <192)1 0 FILES? 4 REC IN USEf 3720 LEFT (of 3724)» OX FULL <14 CYL)» 3330f R/W 
LARSLIB 29C HAS BEEN ATTACHED AS 19C. 

Y <190 R/0 


19C HAS BEEN LOGGED IN 
EXECUTION BEGINS... 
CLASY STARTED 

AS Y 

DISK. 



PROPORTION RELATIVE TO 
00-00 
01-00 

02-45 03-55 

TOP 

LEVEL * 

1 .000000 

1 

PROPORTION RELATIVE TO 

TOP 

LEVEL = 

1.000000 

1 

PROPORTION RELATIVE TO 

TOP 

LEVEL = 

0.750047 

2 

PROPORTION RELATIVE TO TOP 
INDEX = 2 SYMBOL * 

LEVEL = 

****** 

1 * 000000 

1 


KL"*** INDEX <KL>* 2 NUMBER- 2 RELFRP M 0.0 0.0 

PROPORTION RELATIVE TO TOP LEVEL = 0.816835 2 

00*00 
01-00 

02-82 03*18 

04-37 05-45 


PROPORTION RELATIVE TO TOP LEVEL =■ 1.000000 1 

0***SEPERATE 1 SUPER f SUBS 0 2 SPFAC 0.18557E 02 

00-00 

02-82 03-18 

04-37 05-45 


PROPORTION RELATIVE TO TOP LEV'LL = 0.830760 2 

PROPORTION RELATIVE TO TOP LEVEL = 0.484035 5 

00-00 

02-82 03-16 

04-34 05-48 

06-24 07-24 


PROPORTION RELATIVE TO TOP LEVEL * 0.852994 

0***SEPERATE 2 SUPER f SUBS 0 4 SPFAC 0.53649E 02 

00-00 

04-29 05-56 03-15 

06-45 07-11 


rt» J SSv AL PAGE ls 

POOR QUALITY 
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■ ■*** ‘fiT litw'iiwl' 


>twpe classwn exec 


CLASS' YN EXEC FI LI 


4 CONTROL OFF 

CF‘ SPOOL D CONT NOH TO BATCH 
31 PUNCH BATCH MACHINE BATJSC 
& PUNCH BATCH ID J8C23S J8C238 PENT L. 

31 PUNCH BATCH OUTPUT HOUSTON HOUSTON 
4PUNCH EXEC*$ 

4 PUNCH CP SPOOL 9 CLOSE 

& PUNCH CP SPOOL 9 TO JSC 23 5 START CL C NOHOLD 
4 PUNCH GET BISK JSC770 191 291 R B PASS AUCO.IN 
4 PUNCH GETBISK JSC235 191 39:1 W A PASS WRTE7 
4PUNCH EXEC CFILE 41 42 43 44 45 
4 PUNCH RELEASE B <IO 
4 PUNCH CP Q T 

4 PUNCH CP SPOOL 9 CLOSE SI OP 
4 PUNCH *♦ 

CP SPOOL PUNCH NOCONT CLOSE 
CP SPOOL PUNCH HOLD TO RSCS 
4 EX I T 


>tvpp dot e>;ec 


DOT EXEC FILE 




tCONTROL OFF 

GETDISK .'3C770 191 350 R B/A PASS AUCOIN 
GETDISK TEMP 3M CLEAR DETACH 
FIL 21 DISK DOT CC A (PLKSIZE 80 RECFM F LRECL 80 PERM 

FIL 22 DISK FILE 'T22F001 D (LRECL 800 BLKSIZE 800 XTENT 2100 PERM 

FIL 3 TERMINAL (PERM 
FIL 6 PRINTER ( PERM 

FIL FT11F001 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F002 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F003 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F004 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F005 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F0CA TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F007 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F008 TAPI (BLKSIZE 12356 RECFM U PERM 

PIL FT11F009 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F010 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F011 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F012 TAPI (BLKSIZE 12356 RECFM U PERM 

e '•tufOIS TAPI (BLKSIZE 12356 RECFM U PERM 

' "AM TAPI (BLKSIZE 12356 RECFM U PERM 

FIL . - i API (BLKSIZE 12356 RECFM U PERM 

FIL FI a IF 016 TAPI (BLKSIZE 12356 RECFM 'J PERM 

FIL FT11F017 TAF’l (BLKSIZE J2356 RECFM U PERM 

FIL FT11F018 TAF'l (BLKSIZE *2356 RECFM U PERM 

FIL FT11F019 TAF’l (BLKSIZE 12356 RECFM U PERM 

FIL FT 1 1F020 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F021 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F022 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F023 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F024 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F025 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT11F026 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT 1 1F027 TAPI (BLKSIZE 12356 RECFM U PERM 

FIL FT12F001 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F002 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F003 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F004 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F005 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F006 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT 1 2F007 TAP2 (BLKSIZE 12356 PECFM U PERM 

FIL FT12F008 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F009 TAF'2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F010 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F01 1 TAF‘2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F012 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F013 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F014 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL F712F015 7AF2 (BLKSIZE 12356 RECFM U PERM 

FIL FT 12F016 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F017 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F018 TAF‘2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F019 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F020 TAF‘2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F021 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F022 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F023 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT 12F024 TAF‘2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F025 TAF’2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F026 TAF‘2 (BLKSIZE 12356 RECFM U PERM 

FIL FT12F027 TAP2 (BLKSIZE 12356 RECFM U PERM 

FIL FT19F001 DISK 81 D0T2 A (LRECL 80 PLKSIZE 80 PERM 

FIL FT19F002 DISK 81 D0T1 A (LRECL 80 PLKSIZE 80 PERM 

LOAD MONTOR BLKCOM (CLEAR START NOMAP 
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SAMPLE DOT CONTROL CARD IVLi 


>TYPE DOT CC A 


♦DOTDATA 

CHANNEL 

DATAFILE 

DOTFIL 


OPTION 


OPTION 


*END 


DOT 2 

N 

DOT 2 

N 

DOT 2 

N 

DOT 2 

N 

DOT 2 

N 

DOT 2 

N 

DOT 2 

N 

DOT 2 

N 

DOT 2 

N 

DOT 2 

N 

DOT 2 

S 

DOT 2 

S 

DOT 2 

S 

DOT 2 

S 

DOT 2 

S 

♦END* 


♦EXIT 



BATA=l>2,3»4»5>6>7f8 
UNIT- 11 »PIl.E-7 
OUTPUT /UNIT -^l 9 » F 1 LE~1 
LAC IE 
PRINT 

1 6 7 8 10 11 13 14 1? 21 25 29 33 35 37 
39 40 41 42 44 45 47 49 50 SI S3 54 55 

56 57 59 61 65 67 73 78 83 84 86 87 88 9.1 

92 94 95 99 101 105 109 111 113 121 122 123 125 

128 129 130 131 132 135 139 141 145 147 149 153 154 

155 156 159 160 161 162 163 165 1 66 167 168 169 170 

171 175 177 179 183 185 189 191 192 193 194 195 196 

197 203 204 207 208 209 20 22 28 32 34 36 38 
60 62 64 68 72 100 104 136 140 142 146 152 172 
174 176 180 182 186 190 

23459 12 15 16 17 18 23 27 31 43 46 48 
52 63 69 71 75 77 79 80 81 82 85 89 90 93 97 103 
107 115 116 117 118 119 120 124 126 127 133 137 143 151 

157 158 164 173 181 187 198 199 200 201 202 205 206 

26 30 58 66 70 76 96 102 106 134 148 150 178 188 


SAMPLE DOT EXECUTION of Seqment 1394 

>def stor 2i»i 

STORAGE = 02048K 

R» T=0.01/0.02 10J33J06 


>dot 1394 

JSC770 191 HAS BEEN ATTACHED AS 350. 

B (350) R/0 

350 HAS BEEN LOGGED IN AS B/A DISK. 

DASD 192 DETACHED 

YOU ALREADY HAD A VIRTUAL DEVICE 192. IT IS BEING DETACHED. 
TEMP 155 HAS BEEN ATTACHED AS 192. (003.00 MEGABYTES) 

192 HAS BEEN LOGGED IN AS D DISK. 

THE FOLLOWING NAMES ARE UNDEFINED! 

WRTBM 

EXECUTION BEGINS. .♦ 
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>tvpe label exec LABEL EXEC FILL 

$ CONTROL OFF 

PETBISK JSC770 191 350 R B/A PASS AUCOIN 
I BEEP - 1 

-GET GETDISK TEHF* 3M CLEAR DETACH 

ft IF IRETCODE EO 0 ftSKXP 7 

I IF ftDEEP LE 7 ftSKIP 3 

I TYPE NO 10 CYL TEMP DISK AVAILABLE 

CP 0 T 

tEXXT 

CP SLEEP 5 MIN 
ftDEEP = l DEEP + 1 
ft GOTO -GET 
GETDISK LARSYS 
0 DISK * 

FILEDEF 3 TERM (PERM 

FILEDEF 6 DISK LADEL6 LISTING D (PERM 


FILEDEF 

8 DISK ftl DOT ] 

. A (LRECL 3060 RECFM V DLKSIZE 

3060 PERM 

FILEDEF 

FT11F001 

TAPI 

(DLKSIZE 

3060 

RECFM 

U 

DEN 

1600 

PERH 

FILEDEF 

FT1IF002 

TAPI 

(DLKSIZE 

3060 

RECFM 

U 

DEN 

1600 

PERM 

FILEDEF 

FT11F003 

TAPI 

(DLKSIZE 

3060 

RECFM 

U 

DEN 

1600 

PERM 

FILEDEF 

FT11F004 

TAPI 

(DLKSIZE 

3060 

RECFM 

U 

DEN 

1600 

PERM 

FILEDEF 

FT11F005 

TAPI 

(DLKSIZE 

3060 

RECFM 

u 

HEN 

1600 

PERM 

FILEDEF 

FT11F006 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 

FILEDEF 

FT11F007 

TAPI 

(DLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F008 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F009 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F010 

TAF1 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F011 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 

FILEDEF 

FT11F012 

TAPI 

(DLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F013 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 

FILEDEF 

FT11F014 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F015 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F016 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F017 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F018 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 

FILEDEF 

FT11F019 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F020 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

DEN 

1600 

PERM 

FILEDEF 

FT11F021 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F022 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FTUF023 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F024 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FTHF025 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F026 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F027 

TAPI 

(DLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F02B 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F029 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

FT11F030 

TAPI 

(BLKSIZE 

3060 

RECFM 

u 

PEN 

1600 

PERM 

FILEDEF 

12 DISK ftl VARM A (LRECL 80 

BLKSIZE 

80 PERM 



► FILEDEF FT19F001 DISK FILE FT19F001 A (LRECL 1860 BLKSIZE i860 PERM 

FILEDEF FT19F002 DISK FILE FT19F002 A < LRECL 1860 BLKSIZE 1860 PERM 
| FILEDEF FT20F001 DISK FILE FT20F001 A (LRECL 320 BLKSIZE 320 PERM 

FILEDEF FT20F002 DISK FILE FT20F002 A (LRECL 320 DLKSIZE 320 PERM 
FILEDEF FT21F001 DISK LADEL CC A (LRECL 80 DLKSIZE 80 PERM 
FILEDEF FT21F002 DISK FILE FT21F002 A (LRECL 80 DLKSIZE 80 PERM 
FILEDEF FT22FOOI DISK FILE FT22F001 A (LRECL BOO DLKSIZE 000 XTENT 2100 PERM 
t FIL FT23F001 DISK ftl 5 A4 (LRECL 120 DLKSIZE 130 RECFM VS PERM 

r* FIL 24 DISK FILE FT24F001 D A (LRECL 1200 DLKSIZE 1250 RECFM VS PERM 

FILEDEF 25 DISK FILE FT25F001 D (LRECL 132 DLKSIZE 132 PERM 

FILEDEF 26 DISK FILE FT26F001 D (LRECL 132 DLKSIZE 132 PERM 

FILEDEF 27 DISK FILE FT27F001 D (LRECL 132 DLKSIZE 132 PERM 

FILEDEF 28 DISK FILE FT28F001 D (LRECL 132 DLKSIZE 132 PERM 

FILEDEF 29 DISK FILE FT29K001 D (LRECL 132 DLKSIZE 132 PERM 

FIL FT31F001 DISK FILE FT31F001 D (LRECL 3060 BLKSIZE 3060 RECFM U PERM 
GLOBAL TXTLIB FORTRAN CMSLIB 
* LOAD IMAGE BLKCLA (CLEAR START NOHAP 

1 LOAD MAXLAPEL BLKCLA (CLEAR START NOMAP 

\\ r CP REMOTE E TO HOUSTON 

rj CP SPOOL PRINTER HOLD NOCONT 

1 CP TAG OUE DEV E 

L * * D (ALL 

, 1 PRINT LABEL6 LISTING D 

ft IF ft INDEX LT 2 I SKIP 8 
1 M ft IF ft2 EO 1 1SKIP 2 

ftSP « ft2 - 1 
TAPE FSF ftSP (TAP2 

FIL INMOVE DISK FILE FT31F001 D (LRECL 3060 BLKSIZE 3060 RECFM U PERM 

FIL OUT MOVE TAP2 ( BLKSIZE 3060 RECFM U BEN 1600 PERM 

MOVEFILE 

TAPE WTM 2 (TAP2 
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>TYPE LABEL CC A 


SAMPLE LABEL CONTROL CARD FI I,? 


CHANNELS 

MAP 

DATA 

SEGMENT 

PROC 

♦END* 

♦END* 


I»2r3,4,5r6>7,8 

PIXEL 

FILE = 8 

1 394 

07 f 02 r 80 


(l»l>> (1>1)» <196rl>? (196*117) * (1*117) 


SAMPLE EXECUTION OF IMAGE AND LABEL 


FOR SEGMENT 


Mabel 1394 

JSC 770 191 HAS BEEN ATTACHED AS 350* 

7 350' REPLACES ' B (350) ' 

B (350) R/0 

350 HAS BEEN LOGGED IN AS B/A DISK ♦ 

DASD 192 DETACHED 

YOU ALREADY HAD A VIRTUAL DEVICE 192, IT IS BEING DETACHED, 
TEMP 155 HAS BEEN ATTACHED AS 192, (003,00 MEGABYTES) 

192 HAS BEEN LOGGED IN AS D DISK. 

LARSLIB 29C HAS BEEN ATTACHED AS 19C, 

Y (190 R/0 

19C HAS BEEN LOGGED IN AS Y DISK, 

EXECUTION BEGINS,,. 

IMAGE DATA EXTRACTED FROM TAPE * LAST LINE = 117 


EXECUTION BEGINS, ♦. 


3.1 HARDWARE DESCRIPTION 


The CLASSY clustering program is operational on the IBM 3031. 
system at Purdue. 

3.2 EXEC FILES 

3.2.1 CLASSY EXEC FILES 

CLASSY is called by using one of the following EXEC files: 

CLASSY (segment number) , (input tape number) , (input file 
number) , (output tape number) , (output file number) 
for interactive runs. 

CLASSYN (segment number) , (input tape number) , (input file 
number) , (output tape number) , (output file number) 
for batch runs. 

Where segment number = name of control card file with type CC 

input tape number = number of input tape 

input file number = number of input file 

output tape number = number of output tape 
output file number = number of output file 

These EXEC files assign (segment number) CC to be the control 
cara input file which specifies the program options. The output 
files are (segment number) 1, ..., (segment number) (last 
iteration) and FT16F001. 

3.2.2 DOT EXEC FILE 

The pixel data files are extracted by executing the following 
EXEC file to call the DOT DATA processor. 

DOT (segment number) . 
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Where segment number = the name of the output files, (segment 
number) D0T1 and (segment number) D0T2. The name of the control 
card file is DOT CC. 

***** Warning: The tape must be mounted and positioned prior 

to running the DOT EXEC. 

3.2.3 IMAGE AND MAXLABEL EXEC FILE 

The post-processor MAXLABEL and its associated program IMAGE 
are called by the following EXEC file: 

LABEL (segment number) 

where segment number is the number of the segment. 

This EXEC file assigns LABEL CC to be the control card input 
file which specifies the program options. The report output is 
sent to the line printer and the labelled cluster map is placed 
on the D disk as FILE FT31F001. 

***** Warning: The tape must be mounted prior to 

to running the LABEL EXEC. 


3.3 CONTROL CARDS 


3.3.1 CLASSY CONTROL CARDS 

The following control cards are input to the modified CLASSY 
program, and are analysed by SETUPS). In all cards the "keyword” 
begins in card column 1 and any parameters on the card are in 
card columns 11 through 72. Numbers in a series are separated 
by commas; blanks are optional. 

1. "CHANNEL" CARD 

EXAMPLE: CHANNEL 1, 5, 9, 13 

The "CHANNEL" card specifies the channel numbers to be used in 
clustering the multi-channel data vectors. The maximum number 
of channels allowed is sixteen. 

2. "NPTS" CARD 

EXAMPLE : NPTS 2 

This card specifies the number of pixels to skip between 
selected pixels. 

3 . "DATA" CARD 

EXAMPLE: DATA FILE = 2 

This card specifies the input file number. The default value 
is 1. 

4. "MAP" CARD 

EXAMPLE: MAP 1, 3, 5 

This card specifies the iteration numbers for which intermediate 
pixel maps are to be drawn. 

5. "LINES" CARD 

EXAMPLE: LINES 1-10, 25-34 

This card specifies the lines to be mapped on the intermediate 
pixel maps. 


6. "SEGMENT" CARD 

EXAMPLE: SEGMENT 1234 

This optional card is used to specify the segment number used 
in the header of FILE FT31F001. The Accuracy Assessment system 
of programs requires this data in the header. 

7. "PROCESS" CARD 

EXAMPLE: PROCESS 2,11,80 

This optional card is used to specify the date used in the 
header of FILE FT31F001. The Accuracy Assessment system of 
programs requires this data in the header. 

8. "*END*" CARD 

This mandatory card specifies the end of the control cards. 


3.1.2 LABEL CONTROL CARDS 

The following control cards are input to the post-processor 
MAXLABEL and its associated program IMAGE. The cards are 
analysed by the subroutine SETUPM. In all cards, the "keyword" 
begins in card column 1 and any parameters are entered from 
card columns 11 through 72. Numbers in a series are separated 
by commas? blanks are optional. 

1. "CHANNEL" CARD 

EXAMPLE: CHANNEL 1, 5, 9, 13 

The "CHANNEL" card specifies the numbers of the channels to be 
used. The maximum number of channels allowed is sixteen. 

2. "DATA" CARD 

EXAMPLE: DATA FILE =2 

This card is used to specify the input tape file number. The 
default value is 1. 
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3. "MAPOPT" CARD 

EXAMPLES: MAPOPT PIXEL 

MAPOPT CLUSTER 

This card is used to specify either a Pixel Map or a Cluster Map. 

4. "SEGMENT" CARD 

EXAMPLE: SEGMENT 1234 

This optional card is used to specify the segment number used in 
the header of FILE FT31F001. The Accuracy Assessment system of 
programs requires this data in the header. 

5. "PROCESS" 

EXAMPLE: PROCESS 2,1,80 

This optional card is used to specify the date used in the header 
of FILEFT31F001. The Accuracy Assessment system of programs 
requires this data in the header. 

6. "*END*" CARD 

This manditory card specifies the end of the control cards . 

3.3.3 FIELD DEFINITION CARDS FOR CLASSY AND MAXLABEL 
1. " FLDNAM" CARD 

The field definition cards (s) delineate the area on the image 
data tape to be used by CLASSY and MAXLABEL in terms of pixel 
coordinates (sample, line) for each vertex of the "field" up to 
a maximum of 10 vertices for a given field. An alphanumeric 
field identification may be supplied in card columns 1-6 but is 
not required. Coordinate paris are in card column 11 through 72 
and are enclosed in parentheses with the pairs separated by commas. 

The first pair given for a field must be the incrementation 
desired in the lines and pixels to be read from the input image 
tape. I.E., "(2,3)" would indicate every second pixel on each 
line and every third line to be read. 


The second and succeeding coordinate pairs are the (sample, line) 
coordinates of the vertices of the field. A continuation of 
coordinate pairs on the next card is indicated by an asterisk 

Up to 10 coordinate pairs (vertices) are accepted for one 

field. 

EXAMPLE: FLDNAM (2,3) , (2,1), (196,3) , (100,50) * 

(196,100) , (1,100) r (20,30) 

The result of the above cards is that a rectangular area is read 
from the image data which bounds the given irregularly shaped 
field defined above. The coordinates of the rectangular area are 
(1,1) , (196,1) , (196,100) , (1,100). 

The coordinates for the rectangular area are self-determined 
by the tape (file) reading program. Internally, the actual 
field coordinates which were input are used to extract only the 
pixels that are within the actual field defined. 

The input field vertices must be defined on the card(s) in 
clockwise order . 

2. "$END*" CARD 

This mandatory card specifies the end of the field definition 
cards . 
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3. 4 SOFTWARE DESCRIPTION 


3.4.1 CLASSY SUBROUTINES MODIFIED 

3. 4. 1.1 Software Component No. 1 (CLINIT) 

Linkage 

CLINIT is called from CLASSY. 

Interface 

Interface is accomplished through a calling argument and the 
following common blocks: 

/CLUS/, /MISC/, /STPAR/ , /CLUSTR/, /INITL/ and /MXLL/ . 

Input 

KROT - Index to root link. 

Output 

Common blocks initialized. 

Storage Requirements 
N/A 


3. 4. 1.2 SOFTWARE COMPONENT NO. 2 (SETUP9) 

Purpose 

Read control cards for CLASSY. 

Linkages 

Interface is accomplished through common blocks /INFORM/, 
/CLUSTR/ , /FILE/, /MAP/ /TIMERR/, and /WRTAP/. 


Outputs 




NOFSKP 

/FILE/ 

Incremental number of files 

to skip 

NOFEAT 

/INFORM/ 

Number of channels 


FETVEC 

/INFORM/ 

Vertices array 


NPTS 

/CLUSTR/ 

Number of pixels to skip between pixels 

NOCYCL 

/CLUSTR/ 

Number of iterations 


MAP 

/MAP/ 

Array of iterations to map 


MAPCT 

/MAP/ 

Count of maps 


LINES 

/MAP/ 

Matrix of lines to map for 

each iteration 

LINECT 

/MAP/ 

Count of lines to map 


TIMEMX 

/TIMERR/ 

Maximum time 


VAR3L 

/WRTAP/ 

Segment number and date 



Storage Requirements 
Not applicable. 


Description 


SETUP9 reads the control cards in CLASSY CC and saves the 
segment, date, tape and file numbers, map description, lines 
description, maximum time and number of iterations. 


Flowchart 
Not applicable. 

Listing 

See Appendix A for program. 
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3.4.1. 3 SOFTWARE COMPONENT NO. 3 (STATIS) 


Purpose 

STATIS is the control subroutines for the computation and 
reporting part of CLASSY. 

Linkages 

STATIS is called by MULTI, STATIS calls DISC, CLASY2, C0RJ3CT, 
DOTSQ, VPV, VMTV, MPVS, ADJUST, CLDUMP and EXP. 

Interface 

Interface is accomplished through calling arguments and the 
following common blocks: /CLUS/, /MISC/, /STPAR/ , /CLUSTER/, /RAND/, 
and /MXLL/ . 

I nputs 

Temporary scrambled pixel file created by READTP. 

Outputs 

STATIS calls CLDUMP to output the data. 

Storage Requirement 
Not applicable 

Description 

STATIS was modified for the Maximum Likelihood Program to cause 
data to be calculated and saved at the end of each iteration 
through the complete data set. The subroutine additions (1) 
call CALRPR to calculate the relative proportions, and (2-) set 
MXLL in common block /MXLL/ to 1 to cause CLPR to write a file 
record for each active cluster. 
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3.4. 1,4 SOFTWARE COMPONENT NO. 4 (CLPR) 

Purpose 

Report intermediate data to a print file; save data for post- 
processor MAXLABEL. 

Linkages 

CLPR is called from CLDUMP, ADJUST, SEPER and JOIN. CLPM calls 
LOCK, MORSTR, FREE, MORSTR and MINV. 


Interface 

Interface is accomplished through a calling sequence and the 
following common blocks /CLUS/, /MISC/, /STPAR/ and /MXLL/. 


Inputs 


KL - Index to LINK 


IN - .Not used 
SUM - Data to be 
SKEW - Data to be 
KURT - Data to be 
RELP ~ Data to be 
VOLRT - Data to be 
DCON - Data to be 
MEANS - Data to be 
COVAR - Data to be 


displayed 
displayed 
displayed 
saved on file 
saved on file 
saved on file 
saved on file 
saved on file 


Outputs 


Report on report file. 


Storage Requirement 
Not applicable. 
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Description 

CLPR reports to all of the mathematical variables on a report 
file. The variables RELP, VOLRT and DCON and the MEANS array 
and COVAR matrix are saved on the file for the post-processor 
MAXLABEL. 

Flowchart 
Not applicable. 

Listing 

See Appendix A for program. 
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3.4.2 CLASSY SUBROUTINE ADDED 

3. 4. 2.1 SOFTWARE COMPONENT NO. 1 (CALRPR) 

Purpose 

Calculate the relative probability of each pixel belonging to 
each class. 

Linkages 

CALRPR is called by STATIS CALRPR calls ISPLIT. 

Interface 

Interface is accomplished through the /CLUS/, /MISC/, /STPAR/ 
and /MXLL/ common blocks. 


Inputs 





KROTIN 

Calling sequence 

Index 

to 

root link 

INDEX 

/CLUST/ 

index 

to 

cluster data 

LSUBS 

/CLUST/ 

chain 

of 

sub-cluster 

PROP 

/CLUST/ 




PRIRCM 

/CLUST/ 




Outputs 






RELPRP /MXLL/ Relative proportions array 

Storage Requirement 
Not applicable. 

Description 

CALRPR determines the relative proportion for each link from 
the root link by dividing the PROP of the current link by the 
product of PRIRCM for the previous link and the. relative proportion 
for the previous link. 
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Flowchart 
Not applicable. 

Listing 


See Appendix A for program 


3.4.3 IMAGE PROCESSOR AND SUBROUTINES 


3.4. 3.1 SOFTWARE COMPONENT NO. 1 PROCESSOR (IMAGE) 

Purpose 

Create an input file for the MAXLABEL processor from a JSC 
Universal Image Tape. 

Linkages 

Interface is accomplished through common blocks /INFORM/ and 
/CLUSTR/. 

Outputs 

Pixel file for MAXLABEL Processor. 

Storage Requirement 
Not applicable. 

Description 

IMAGE calls SETUP9 to read the LABEL control cards to describe the 
the pixel data needed for the MAXLABEL processor. The subroutine 
KREDTP calls the LARSYS subroutines to read the JSC Universal 
Image Tape and moves the pixel data to the pixel file for the 
MAXLABEL processor. 

Flowchart 
Not applicable. 

Listing 

See Appendix A for program. 
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3.4. 3. 2 SOFTWARE COMPONENT NO. 2 (SETUPM) 

Purpose 

Read control cards for IMAGE and MAXLABEL. 

Linkages 

Interface is accomplished through common blocks /INFORM/. 
/CLUSTR/ and /WRTAP/. 


Outputs 




NOFSKP 

/FILE/ 

Incremental number 

of files to skip 

NOFEAT 

/INFORM/ 

Number of channels 


FETVEC 

/INFORM/ 

Vertices array 


NPTS 

/CLUSTER/ 

Number of pixels to skip between pixels 

VARBL 

/WRTAP/ 

Segment number and 

date 


Storage Requirements 
Not applicable. 

Description 

SETUPM reads the control cards in LABEL CC and saves the segment, 
date, channel information and map option. 

Flowchart 
Not applicable. 

Listing 

See Appendix A for program. 


3.4. 3.3 SOFTWARE COMPONENT NO. 3 (KREDTP) 


Purpose 

Read lines of data from image tape. 

Linkages 

KREDTP is called the post-processor IMAGE. This subroutine calls 
TAPHDR, FLDINT, LINERD and FDLINT . 

Interface 

Interface is accomplished through common blocks /INFORM/ , 

/CLUSTR/ and /FILE/. 

Inputs 

NOFSKP /FILE/ Incremental number of files to skip 
NPTS /INFORM/ Number of pixels to skip between pixels 

NOFEAT /INFORM/ Number of channels 

FETVEC /INFORM/ Vertices array 

Outputs 

FILE of pixels as described by control cards. 

Storage Requirement 
Not applicable. 

Description 

KREDTP is calls TAPHDR to read the tape header, LAREAD to read 
the field and vertices information, FLDINT to position tape for 
this field, LINERD to read lines from the universal format tape. 

Flowchart 
Not applicable. 
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Listing 

See Appendix A for program. 

Error Message 

FIELD DEFINATION INFORMATION EXCEEDS 2000 WORDS. END-OF-TAPE 
REACHED BEFORE END OF FIELD. 

Data Mesaages 


Vertices listed 


3.4.4 MAXLABEL PROCESSOR AND SUBROUTINES 


3. 4. 4. I SOFTWARE COMPONENT NO. 1 (MAXLABEL) 

Purpose 

Maxlabel estimates the probability of observing a particular 

labelled class given that a particular cluster has been observed 

(BETA) . These estimates are used to produce estimates of the 
• ■ 

proportion of each labelled class present in the scene. The 
proportion estimates are computed in two different ways; The 
first assumes that clusters are to be "bias corrected" using the 
BETA'S. The second assumes that clusters are labelled using the 
BETA s. 

Linkages 

MAXLABEL calls SETUPM, ALLPXI , BAPLS , PRTBAP, PRTAP , PRTAP2 and 
WRTLNS . 

Interface 

Interface is accomplished through calling sequences, blank common 
and /FILE/ common blocks. 

Inputs . 

Control caird file LABEL CC 

Image input file 

CLASSY cluster statisical parameters. (segment number) 1, ..., 

(segment number) (last iteration) 


Outputs 

Iteration Report (Estimated BETA'S, cluster labels 

and propotioris of each class present.) 

Pixel Labelled Dot Map 
Cluster Labelled Dot Map 
Unlabelled Dot Map 

Entropy and Posterior Probability for each Cluster Report 
Pixel Labelled Cluster Map 
Cluster Labelled Cluster Map 
Un labelled Cluster Map on Tape 

Storage Requirement 
Not applicable. 

Description 

MAXLABEL (1) reads the LABEL CC file to determine the segment 
number, processing date and data description (2) reads the dot 
data file and (3) reads the CLASSY cluster statistical parameters 
file. The iterative process then begins to calculate the 
estimates. The class and cluster labels are saved for maximum 
estimates. 

Flowchart 
Not applicable. 

Listing 

See Appendix A for program. 
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3. 4. 4. 2 SOFTWARE COMPONENT NO. 2 (SETUPM) 

See 3,4, 3.2 for description of SETUPM. 

3. 4. 4. ,3 SOFTWARE COMPONENT NO. 3 (READCC) 

Purpose 

Read CLASSY cluster statistic parameters into common block /CLASY/. 
Linkage 

READCC is called from MAXLABEL. 


Interface 

Interface is accomplished through blank common and /CLASY/ blocks. 


Inputs 

File from CLASSY. (Segment number) (Last iteration number) . 


Outputs 


NOCC 

RLPRP 

CCLRT 

CDCON 


/blank common/ 
/blank common/ 
/blank common/ 
/blank common/ 


Number of CLASSY clusters 
Relative proportion array 
Normalization factor elements 
Normalization factor elements 


Storage Requirements 


Not applicable. 

Description 

READCC reads the data from CLASSY describing the clusters. 


Flowchart 
Not applicable. 

Listing 

See Appendix B for program. 

3-33 


mi Tuiirfuinra nr 




3. 4. 4. 4 SOFTWARE COMPONENT NO. 4 (ALLPXI) 

Purpose 

Calculate the probability of X (dot data pixel) given I (Clasy 
Cluster) for each X in the dot data file. 


Linkages 

ALLPXI is called by the main program MAXLABEL. ALLPRI calls 
DOTSQK and GETCC . 

Interface 

Interface is accomplished through blank common. 


Inputs 

NOCC 

/blank 

common/ 

ITOTDT 

/blank 

common/ 

MQ 

/blank 

common/ 

IDOTS 

/blank 

common/ 

CMEANS 

/blank 

common/ 


Number of CLASSY Cluster 
Total number of dots 
Number of channels 
Dot data pixels 
CLASSY cluster mean matrix 


Outputs 


PXI /blank common/ Probability of pixel belonging to CLASSY 

cluster. 


Storage Requirement 
Not applicable. 

Description 

The array PXI is created by subtracting the CLASSY cluster mean 
value from the dot data pixel value for each channel . The product 
of this array and the covariance matrix for the CLASSY cluster 
is created. The probability is calculated as P(X.I) = E 
( -PRODUCT/2.) *E (-DC0N/2 . ) /CVOLRT 
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Where DCON and CVOLRT are the normalization factor elements of 
the CLASSY cluster. 

Flowchart 
Not applicable. 

Listing 

See Appendix A for program. 


\ 
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3. 4. 4. 5 SOFTWARE COMPONENT NO. 5 (BAPLS) 

Purpose 

Sum products in a class for all dot-cluster combinations. 
Linkages 

BAPLS is called by MAXLABEL . 

Interface 

Interface is accomplished through blank common. 


Inputs 




NOCC 

/blank 

common/ 

Number of CLASSY cluster 

NOCAT 

/blank 

common/ 

Number of classes 

ITOTDT 

/blank 

common/ 

Total number of dots 

IDOTS 

/blank 

common/ 

Dot data pixels 

BETA 

/blank 

common/ 

Prior estiamte 

RLPRP 

/blank 

common/ 

Relative porportion for cluster 

PX 

/blank 

common/ 

Probability of pixel belonging to 


cluster 


Outputs 

SLK /blank common/ Sum of products by class and cluster 
SK /blank common/ Sum of products by cluster 

Storage Requirement 
Not applicable. 

Description 

For each dot the following calculations are made: 

(1) The class of the dot is determined 

(2) The products of BETA (class, cluster) *RLPRP (cluster) 
*PX (cluster, dot) are summed. 


Flowchart 


Not applicable. 

Listing 

See Appendix A for program. 
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3. 4. 4. 6 SOFTWARE COMPONENT NO. 6 (PRTELB) 

Purpose 

Construct and print estimate for each labelled class. 
Linkages 

PRTELB is called by the main program MAXLABEL. 

Interface 

Interface is accomplished through blank common. 

Inputs 

BETA /blank common/ Prior estimate 

RLRP /blank common/ Relative proportion 

Outputs 

Report on report file. 

Storage Requirement 
Not applicable. 

Description 

PRTELP creates the sum of the products BETA*RLPRP for each 
class-category combination. 

Flowchart 
Not applicable. 

Listing 

See Appendix A for program. 
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3. 4. 4. 7 SOFTWARE COMPONENT NO. 7 (PRTBAP) 

Purpose 

Construct and print labeled cluster map for each class-cluster. 
Linkages 

FRTBAP is called by the main program MAXLABEL. 

Interface 

Interface is accomplished through the blank common and /MPPXL/ 
common blocks . 

Inputs 


ITOTDT 

/blank common/ 

Total dots 

IDOTS 

/blank common/ 

Dot location 

NOCAT 

/blank common/ 

Number of categories 

BETA 

/blank common/ 

Prior estimate 

RLPRP 

/blank common/ 

Relative proportion 

PX 

/blank common/ 

Pixel value 

LABELS 

/blank common/ 

Category labels 


Outputs 

MPXLA /MPPXL/ Pixel labels. 

Storage Requirement 
Not applicable . - 

Description 

The maximum product for each prior estimate * relative proportion 
* pixel value is determined. A matrix is constructed of the 
class labels corresponding to the maximum products . This matrix 
is printed. 
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Flowchart 


Not applicable. 

Listing 

See Appendix A for program. 
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3. 4. 4. 8 SOFTWARE COMPONENT NO. 8 (PRTAP) 

Purpose 

Construct and print labeled cluster map for each cluster. 


Linkages 

PRTAP is called by the main program MAXLABEL. 


Interface 

Interface is accomplished through blank common and /MPPXL/ 
common blocks . 


Inputs 

ITOTDT 

RLRPR 

PX 

NOCC 
LB LOST 


/blank common/ Total number of dots 
/blank common/ Relative proportion for class 
/blank common/ Pixel value for class 
/blank common/ Number of CLASSY cluster 
/blank common/ Labels for CLASSY clusters 


Outputs 

MPXLA /MPPXL/ Class labels 


Storage Requirement 
Not applicable. 

Description 

The maximum product for each relative proportion * pixel value is 
determined. A matrix is constructed of the class labels 
corresponding to the maximum products. This matrix is printed. 
The entropy for each class-CLASSY cluster and the probability 
for each pixel-CLASSY cluster are calculated and printed. 
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Flowchart 


Not applicable. 

List ing 

See Appendix A for program. 
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3. 4. 4. 9 SOFTWARE COMPONENT NO. 9 (PRTAP2) 
Purpose 

Print labelled cluser map for each cluster. 


Linkages 

PRTAP2 is called by the main program MAXLABEL. 


Interface 

Interface is accomplished through blank common and /MXPPXL/ 
common blocks. 


inputs 


ITOTDT 

NOCC 

RLPRP 


LBLCSS 


/blank common/ Total number of dots 
/blank common/ Number of CLASSY clusters 
/blank common/ Relative proportion for CLAS3Y cluster 
/blank common/ Pixel value for CLASSY cluster 
/blank common/ Labels for CLASSY cluster categories 


Outputs 

Line written to report file of CLASSY cluster categories 
corresponding to pixel. 

Storage Requirement 
Not applicable. 

Description 

The maximum product for each relative proportion * pixel value 
is determined. A matrix is constructed of the cluster labels 
corresponding to the maximum products. This matrix is printed. 
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Flowchart 


Not applicable. 

Listing 

See Appendix A for program. 


3.4.4.10 SOFTWARE COMPONENT NO. 10 (WRTLNS) 

Purpose 

(1) Read radiance values for lines of original pixel data, 

(2) Print heading for reports, and 

(3) Call subroutines to calculate report data. 

Interface 

WRTLNS is called by MAXLABEL. WRTLNS calls WRTHED, PXILN, LNBAP, 
WRTLN, LNAP, LNAP2 PAGE. 

Inputs 

Temporary file on unit 24. 

Outputs 

Pixel Labelled Cluster Map 
Cluster Labelled Cluster Map 
Unlabelled Cluster Map 

Storage Requirement 
Not applicable. 

Description 

WRTLNS reads a description of the lines from the temporary file. 
Report headings are written to files to be saved for each of 
the reports. The lines of radiance values are then read, the 
report calculations made and the report lines written to 
temporary report files. The temporary report files are then 
sent to the printer. 
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Flowchart 


Not applicable. 
Listing 


See Appendix A for program 


3.4.4.11 SOFTWARE COMPONENT NO. 11 ( LNBAP ) 


Purpose 

Print labelled cluster map for each class-cluster. 

Pixel cluster = Max (Beta * alpha * P(X.I)). 

Linkages 

Interface is accomplished through /PXLLN/ , /MPPXL/ and blank 
common blocks. 

inputs 

IBEGIN /PXLLN/ 

IEND /PXLLN/ 

BETA /blank common/ 

RLPRP /blank common/ 

PXLN /PXLLN/ 

Outputs 

Report line written on labelled cluster map file. 

Storage Requirement 
Net applicable. 

Description 

Each pixel in the line of data is processed for every category- 
cluster combination to determine the maximum sum of the class 
products for the category where the product = the estimate * 
the relative proportion * the probability that the pixel belongs 
to the cluster. 


First pixel on line 

Last pixel on line 

Prior estimate 

Relative proportion 

Probability of pixel given cluster 


Flowchart 


Not applicable. 

Listing 

See Appendix A for program. 
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3.4.4.12 SOFTWARE COMPONENT NO. 12 (PXILN) 

Purpose 

Calculate the probability of X (dot data pixel) given I (CLASSY 
cluster) for each X in the line in the image file. 

Linkages 

PXILN is called by WRTLNS . PXILN calls GETCC and DOTSQK. 
Interface 

Interface is accomplished through blank common and /PXLLN/ . 


Inputs 

* LNDOTS 
CMEANS 
MQ 

CVRIN 

/PXLLN/ 

/blank common/ 
/blank common/ 
/blank common/ 

CLASSY mean array 
Number of channels 
CLASSY covariance matrix 

DCON 

/blank common/ 

Normalization factor 

element 

CV'OLRT 

/blank common/ 

Normalization factor 

element 

Outputs 

PXLN /PXLLN/ Probability of X given cluster. 



Storage Requirements 
Not applicable. 

Description 

PXILN makes the pixel data calcualtions for each CLASSY cluster 
by calling GETCC to read the data for the cluster and then 
processing each dot data pixel using that cluster data to 
calculate the probability of that dot data pixel given that 
CLASSY cluster. 
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Flowchart 
Not applicable. 




See Appendix A for program. 




3.4.4.13 SOFTWARE COMPONENT NO. 13 (GETCC) 

Purpose 

Move one set of CLASSY data from arrays to corresponding scalars 


Linkage 

GETCC is called from GETCC. 

Interface 

Interface is accomplished through blank common and /CLASSY/ 
blocks . 


Input 

ICC /blank common/ current class 


Output 


RfeLPRP 

CVOLRT 

CMEANS 

CVRIN 

DCON 


/blank common/ 
/blank common/ 
/blank common/ 
/blank common/ 
/blank common/ 


Relative proportion 
Normalization factor element 
CLASSY cluster mean array 
Covariance matrix 
Normalization factor element 


Storage Requirement 
Not applicable. 


Description 

GETCC moves data for this specified class from the arrays of 
CLASSY data. 


Flowchart 

Not applicable. 

Listing 


See Appendix A for program 
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3.4.4.14 SOFTWARE COMPONENT NO. 14 (DOTSQK) 

Purpose 

Calculate the inner product DIFXMN. DIFXMN relative to the metric 
CVRIN. 

Linkages 

DOTSQK is called from ALLPXI. 

Interface 

Interface is accomplished through calling arguments. 

\ 

Inputs 

MQ Number of channels 

DIFXMN Array of differences between pixel values and CLASSY 
mean. 

CVRIN CLASSY covariance matrix. 

Outputs 

DOTSQK enter product DIFXMN. DIFXMN relative to the metric 
CVRIN. 

Storage Requirement 
Not applicable. 

Description 

DOTSQK calculates the inner product DIFXMN. DIFXMN relative to 
the metric CVRIN. 



3.4.4.15 SOFTWARE COMPONENT NO. 15 (PAGE) 

Purpose 

Read data from temporary file and write report heading and 
column headings. 

Linkages 

PAGE is called by WRTLNS. 

Interface 

Interface is accomplished through a temporary file. 

Inputs 

I UNIT temporary file unit. 

Outputs 

Report heading written on report file. 

Storage Requirement 
Not applicable. 

Description 

PAGE reads the temporary file and writes the data on the report 
file. 

Flowchart 
Not applicable. 

Listing 

See Appendix A for program. 
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4. OPERATION 


CLASSY, DOTDATA, IMAGE and LABEL are operational on the IBM 3031 computer at 
LARS, West Lafayette, Indiana. 

The programs, EXEC files and CC files can be loaded from tape 3956, file 5. 

I 

CLASSY is executed by entering the following commands after signing on the 
computer system. 

OEF STOR 2M 

TAPMOUNT (Tape Number) TAPI RO 1600 
CLASSY or ClASSYN (Segment Number) 

Control Inputs is read from (Segment Number) CC. 

Text output in on the terminal and line printer which are assigned in the EXEC, 
and the One Channel Unlabelled Tape File is on the output tape designated in the 
execute statement. 

DOTDATA is executed by entering the following commands: 

TAPE REW 

DOT (Segment Number) 

Control input is read from DOT CC. 

Text output is on the terminal and line printer, and the Ground Truth Data 
File is on (Segment Number) D0T1 and (Segment Number) D0T2. 

IMAGE and LABEL are executed by entering the following commands: 

TAPE REW 

LABEL (Segment Number) 

Control input is read from LABEL CC. 

Text output is on the terminal and line printer, and the One Channel Labelled 
Tape File is on the tape designated in the execute statement. 


4-1 




CLASSY SETUP9 NXTCHR 

OUMBER 

READTP RINIT 
— TAPHDR 
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— LI NERD 
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— RREAO 
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APPENDIX A 


LISTINGS OF MODIFIED CLASSY SUBPROGRAMS 
LISTINGS OF IMAGE AND MAXLABEL PROGRAMS 
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